THREE TYPES OF SELF-SIMILAR BLOW-UP FOR THE 
FOURTH-ORDER p-LAPLACIAN EQUATION WITH SOURCE: 
VARIATIONAL AND BRANCHING APPROACHES 



V.A. GALAKTIONOV 

Abstract. Self-similar blow-up behaviour for the fourth-order quasilinear p-Laplacian 
equation with source, 

u t = — {\u xx \ n u xx ) xx + \u\ p ' 1 u in M x M + , where n > 0, p > 1, 

is studied. Using variational setting for p = n + 1 and branching techniques for p =/= n + 1, 
finite and countable families of blow-up patterns of the self-similar form 

u s (x,t) = (T-t)-^f(y), where y = x/(T-t)?, (3 = - ^g^ , 

are described by an analytic-numerical approach. Three parameter ranges: p = n + I 
(regional), p > n+1 (single point), and 1 < p < n+ 1 (global blow-up) are studied. This 
blow-up model is motivated by the second-order reaction diffusion counterpart 

u t = (\u x \ n u x ) x +u p ( u >°) 

that was studied in the middle of the 1980s, while first results on blow-up of solutions 
were established by Tsutsumi in 1972. 

This paper is an earlier extended preprint of [2"2"] . 



1. Introduction: classic and recent blow-up reaction-diffusion models 

1.1. Classic second-order model. The nonlinear p-Laplacian operator in Mr, 

(1.1) A p u = V • (| Vw| p ~ 2 Vw), with exponents p>l ( V = grad^,) , 

which serves as a natural extension of the Laplacian 

A = A 2 , i.e., for p = 2, 

enters many classic PDEs of mathematical physics. One of the key mathematical ad- 
vantages of the p-Laplacian (II. ip is that it is nonlinear and at the same time remains 
a monotone operator in the L 2 -metric precisely as the linear Laplacian A does. Oper- 
ators such as (II. ip appear in many works on nonlinear parabolic or elliptic PDEs since 
the 1950s; see various examples, references, and applications in Lions' classic book |39j . 
Gradient-dependent nonlinear operators are typical for filtration, combustion (solid fuels), 
and non- Newtonian (dilatable, pseudo-plastic fluids) liquids theory; see [221 P- 428]. 
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Concerning parabolic PDEs admitting blow-up solutions that compose the main subject 
of the present paper, the p-Laplacian also appeared before other well-known nowadays 
porous medium type nonlinearities (see equation (11.41) below). Namely, it is remarkable 
that the first results on blow-up in quasilinear parabolic equations were obtained by 
Tsutsumi in 1972 [50] for the second- order p-Laplacian equation (pLE-2) with source 
posed in a bounded smooth domain Q C K with the zero Dirichlet boundary condition: 

(1.2) ut = V ■ (|Vu| n W) +u p in tlxl + (u > 0), 

where, in comparison with ( 11. ip . we have renamed the exponents by setting n = p — 2 > 
and write the source term as u p . Concerning the structure of blow-up singularities, vari- 
ous countable and finite families of self-similar blow-up patterns for the one- dimensional 
equation of (II. 2p (and also for the radially symmetric version of (11.21) ). 

(1.3) u t = (\u x \ n u x ) x + u p , 

have been known since the middle of the 1980s; see anc ^ °ther related references 

therein. Surprisingly for the author who initiated the study in [261 12Z] , it turned out that 
(jl.3p generates much wider countable and even uncountable families of self-similar blow-up 
patterns than the porous medium equation with source (PME with source) 

(1.4) u t = (u n+1 ) xx + u p (u>0), 

which was studied by Kurdyumov's Russian School on blow-up and localization since the 
beginning of the 1970s; see history, references, and basic results in [4"71 Ch. 4]. 

It is worth mentioning that the set of blow-up similarity solutions of fll .3D . to say 
nothing about non-radial patterns for (jl.2p . is rather complicated (e.g., contains infinite 
countable and even uncountable families of positive solutions for p > n + 1), so there 
are still some difficult open mathematical problems concerning the structure of blow-up 
singularities for (II. 3p . 

Blow-up results for the p-Laplacian equations with source (I1.3P and (jl.2p together with 
Fujita's pioneering study of the semilinear heat equation (1966) [T8] , 

(1.5) u t = Au + u p (p>l), 

are crucial for modern singularity and blow-up theory of nonlinear evolution PDEs. 

Nowadays, blow-up and other singularity formation phenomena for various classes 
of nonlinear evolution PDEs are rather popular in mathematical literature and appli- 
cations in mechanics and physics. It is well established that blow-up phenomena in 
nonlinear PDEs not only present principal evolution patterns of interest in applica- 
tion, but also can give insight into the deep mathematical nature of nonlinear equa- 
tions under consideration and describe general aspects of various fundamental problems 
of existence-nonexistence, uniqueness-nonuniqueness, optimal regularity classes, and ad- 
missible asymptotics of proper solutions. 

To emphasize that this is not an exaggeration, let us mention that, according to the typ- 
ical tools of the possible and already available analysis and proofs, that the two key open 
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PDE/geometry problems of the twentieth and twenty-first century are directly attributed 
to the area of PDE blow-up research: 

Problem (I): Poincare Conjecture (a closed connected 3D manifold is homeomorphic 
to S 3 ) and the general geometrization problem with Perel'man's recent proof by intro- 
ducing two new monotonicity formulae and others to pass through blow-up singularities 
of Ricci flows with surgery (see [11] for a full account of history, references, and recent 
development); and 

Problem (II): Uniqueness or nonuniqueness (and hence nonexistence or existence of 
local small-scale blow-up singularities) in the 3D Navier-Stokes equations^ (one of the 
Millennium Prize Problems for the Clay Institute; see Fefferman |17j). 

There are several monographs [H H3 HH EES EH [201 US] , which are devoted mainly to 
space-time stricture of blow-up singularities in second-order reaction-diffusion PDEs and 
explain the role of blow-up phenomena in general PDE theory. See also [33] presenting 
various exact solutions and some examples of partial singularity analysis of other classes of 
thin film, nonlinear dispersion, and hyperbolic PDEs. In the monograph [40J, a nonlinear 
capacity approach was shown to be efficient to detect conditions of global nonexistence 
for a variety of nonlinear PDEs and systems of different orders and types. 

The questions of the space-time structure and multiplicity of possible blow-up asymp- 
totics represent problems of higher complexity that need another more involved math- 
ematical treatment, which often and still cannot be fully justified rigorously, so a true 
combination of various approaches including enhanced numerics is in great demand. 

1.2. Fourth-order reaction-diffusion equation. In this paper, we study self-similar 
blow-up for the following quasilinear parabolic fourth-order p-Laplacian equation with 
source (pLE-4 with source): 



u t = A (it) = - (\u xx 


nUxx )xx + 


U 


p^u in RxR + , 



where, as above, n > and p > 1. Here, similar to (11.11) . the fourth-order p-Laplacian 
operator, where we set p = n + 2 > 1 (N = 1 in (11.61) ). 

A Pt2 u = -A(\Au\ p - 2 Au) 
is monotone in the metric of L 2 (R N ). For n = 0, (II. 6p reduces to the semilinear equation 

(1.7) u t = -u xxxx + |u| p ~V 

which describes single point blow-up only for all p > 1 and is already known to admit 
various similarity and other blow-up solutions, [8]. Moreover, it is curious that we have 
found quite fruitful to use the analogy with the linear bi-harmonic equation 

(1.8) u t = -u xxxx + u in R x R + , 

which is obtained from (11.61) by both limits n — > and p — > 1. A simple countable 
subset of exponential patterns for (ll.8p is easy to describe on the basis of spectral theory 



See |24j as a most recent survey, where connections with reaction-diffusion theory are discussed. 
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presented in Section 13.31 Eventually, we will detect certain traces of such countable sets 
(the so-called p-branches) of similarity solutions in the nonlinear problem (11 .6p . 

Being involved in the mathematical study of blow-up for the PME with source (11.41) 
from the middle of 1970s and for the pLE-2 with source fll.2jl from the 1980s, the author 
must admit that the study of blow-up patterns for the proposed pLE-4 with source (11.61) 
was quite a challenge and the author did not expect that the necessary mathematics 
should be so dramatically changed to cover approximately the same concepts developed 
twenty or even thirty years earlier. Recall that in (11.61) we just increase by two the order of 
the diffusion operator in comparison with the standard model (II. 3p . However, this makes 
almost all mathematical tools applied before very successfully to (ll.3p almost nonexistent. 

Thus, we consider for (II. 6p the Cauchy problem with given bounded compactly sup- 
ported data 

(1.9) u(x,0) = uo(x) £ Co(M). 

Since the operator A in (II. 6p is potential in the metric of L 2 and the p-Laplacian is also 
a monotone operator there, local existence and uniqueness of a unique weak (continuous) 
solution, which is defined in the standard manner, are not principal issues and follow from 
classic theory of monotone operators; see Lions [39], Ch. 2]. Finite propagation phenomena 
for the PDE (II .6p are proved by energy estimates via Saint- Venant's principle; see |48j . 
references therein, and a survey in [32J. Therefore, there exists the unique local solution 
of the Cauchy problem (II .6p . (11.91) . which is a compactly supported function u(x,t) that 
can blow up in finite time in the sense that 

(1.10) sup^gjg \u(x, t)\ — > +oo as t — > T~ < oo. 

Existence of blow-up in such higher-order quasilinear parabolic equations is a reasonably 
well-understood phenomenon; see references and approaches in [T3"} [T9~j 128] and Mitidieri- 
Pohozaev [ID]. For instance (see [T3] and references therein), it is known that, for the 
equation similar to (II. 6p with the absolute value in the source-term, 

(1.11) u t = -(\u xx \ n u xx ) we + H P m KxR + , 

all nontrivial solutions with data having positive first Fourier coefficient, 

/ uq{x) dx > 0, 

E 

blow-up in finite time in the subcritical Fujita range 

n + l<|)<pi) = n + l + ^| JV=1 = 3n + 5 (n > 0), 
as well as, most probably, in the critical case p = Po , which needs additional study. 

1.3. Layout of the paper: three types of blow-up. In Section El we describe some 
local and rather delicate oscillatory properties of travelling wave solutions near finite 
interfaces. This is the first time, where we face difficult and still non fully justified 
mathematics concerning higher-order degenerate p-Laplacians. Section [3] is devoted to 
the setting of blow-up self-similar solutions and some preliminaries concerning the linear 
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operator with n = (even this issue is not that straightforward and demands essentially 
non self-adjoint theory). 

Further principal difficulties and important mathematical problems in the study of 
such blow-up solutions concern the description and classification of possible types (the 
structure, stability and multiplicity) of blow-up patterns occurring in finite time. Later 
on, we study three classes of similarity blow-up solutions of (jl.6p in the ranges: 

(i) Section HJ p = n + 1, regional blow-up, so the infinite limit (jl.lOp occurs on a 
bounded x-interval; 

(ii) Section [5j p > n+ 1, single point blow-up, so ( jl.lOp happens at a single point, say, 
at x = 0, and then u(x, T~) is bounded for any x^O; and 

(iii) Section O p G (l,n+ 1), global blow-up, and (jl.lOp happens for any x G R (and 
possibly uniformly on any bounded x-interval). 

A similar classification and various single point blow-up patterns of the so-called P-, 
Q-, R-, and S-type for the second-order counterpart (ll.3p have been known since 1980s; 
see [27] and more references and results in [7J. Actually, we show that some concepts of 
the methodology developed in [27J [7J for (11. 3p also apply to the fourth-order reaction- 
diffusion equation (II. 6p . but indeed demand a different and more difficult mathematics. 
Several problems remain open still. It turns out that, in general, the PDE (11.61) admits 
more complicated sets of similarity patterns than the fourth- order porous medium equation 
(PME-4) with source [23] . 

(1.12) u t = -(\u\ n u) xxxx + \u\ p ~ l u. 

The general scheme of blow-up study via variational and branching approaches applies to 
higher-order p-Laplacian PDEs such as the pLE-6 with source (or any 2mth-order one) 

u t = (\u xxx \ n u xxx ) xxx + \ur l u (or ut = (-l) m+1 D™(\D™u\ n D™u) + {u^u). 

1.4. On some other higher-order PDEs with blow-up, extinction, and finite 
interfaces. Blow-up in parabolic PDEs with higher-order diffusion becomes much more 
difficult than for second-order reaction-diffusion equations. Even simpler PDEs such as 
the extended Frank- Kamenetskii equation in one dimension 

(1.13) Ut — u xxxx + e or Ut — u xxxxxx 4- e , 
and their counterparts with power nonlinearities 

(1.14) u t = -u xxxx + |u| p_1 u or u t = u xxxxxx + \u\ p ~ x u, 

revealed several principally new asymptotic blow-up properties demanding novel math- 
ematical approaches; see details in [HI \W\ . Similar difficulties occur for the Semenov- 
Rayleigh-Benard problem with the leading operator of the form 

(1-15) u t = -u xxxx + (3[{u x f} x + e u (/3>0); 

see [35]. The mathematical difficulties in understanding the ODE and PDE blow-up 
patterns increase dramatically with the order of differential diffusion-like operators in the 
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equations. Interesting regional blow-up and oscillatory properties [9] are exhibited by a 
semilinear diffusion equation with "almost linear" logarithmic source term 

(1.16) u t = -u xxxx + u In 4 \u\. 

The above models are semilinear and do not admit blow-up patterns with finite interfaces. 

Concerning quasilinear higher-order PDEs, the interface and blow-up phenomena are 
natural and most well-known for the degenerate unstable thin film equations (TFEs) with 
lower-order terms such as 

(1.17) u t = -{\u\ n u xxx ) x - {\u\ p ~ l u) xx , ut = -V •(M n VAu)±A|u| p - 1 u, 

where n > and p > 1. Equations of this form are known to admit non-negative solutions 
constructed by special sufficiently "singular" parabolic approximations of nonlinear coeffi- 
cients that lead to free-boundary problems. This direction was initiated by the pioneering 
paper [!3J and was continued by many researchers; we refer to [381 E2] and the references 
therein. Blow-up similarity solutions of the fourth-order TFE (11.171) with the unstable 
sign "— " 

(1-18) ut = ~{u n u xxx ) x - (u p ) xx (u > 0), 

have been also well studied and understood; see [5], El [151 SSI E2] , where further references 
on the mathematical properties of the models can be found. Countable sets of blow-up 
patterns for this TFE were described in [T5] . 

Interface and finite-time extinction behaviour, which is described by various similarity 
patterns, occur for other reaction-absorption PDEs such as 

(1.19) U t = -u xxxx - |w| p-1 m 
in the singular parameter range 

(1.20) p G (— |, 1), =>- is not Lipschitz continuous at u = 0, 

so that is not Lipschitz continuous at u = 0; see [21J and references therein. 

We have used a simply looking quasilinear model such as ( \l.6\\ to demonstrate various 
new aspects of higher-order reaction-diffusion blow-up phenomena. The mathematics 
then becomes more difficult than for the second-order PDEs in (11.3j) . where the Maximum 
Principle reveals its full capacity. We do not expect straightforward rigorous justifications 
of several our conclusions and results, and state key open problems when necessary. 

2. Local asymptotic properties of solutions near interfaces 

Here, we describe generic oscillatory behaviour of solutions of (11. 6p close to finite in- 
terfaces. 
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2.1. Local properties of travelling waves: oscillatory profiles for A < 0. We use 

simple TW solutions, 

(2.1) u(x,t) = f(y), y = x-Xt, 

to check generic propagation properties for reaction-diffusion equations involved. In a 
wide class of ID second-order reaction-diffusion parabolic PDEs, the TWs rigorously 
describe the behaviour of finite interfaces for general classes of solutions; see [20l Ch. 7] 
and references therein. 

We use this approach for the fourth-order PDE (11.61) . The ODE for / takes the form 

(2.2) -xf=-(\n n n n +\fr 1 f- 

By a local analysis near the singular point {/ = 0, /' = 0}, it is not difficult to show 
that the higher-order term \f\ p ~ 1 f on the right-hand side is negligible. Therefore, near 
interfaces, assuming that these are propagating, we can consider the simpler equation 

(2.3) (\f"\ n f")' = -f for y>0, /(0) = 0, 

which is obtained on integration once. Here we set A = —1 for propagating waves, 
by scaling. We need to describe its oscillatory solution of changing sign, with zeros 
concentrating at the given interface point y = + . Oscillatory properties of solutions 
are a common feature of related higher-order degenerate ODEs; see pioneering paper by 
Bernis-McLeod jl] for similar fourth-order ODEs. 

It follows from the scaling invariance of (I2.3P that there exist solutions of the form 

(2.4) f(y)=y fl (p(s), s — lay, where \x = > 2 for n > 0, 

where <f{s) is called the oscillatory component of the given solution. Substituting (12.41) 
into (12. 3p yields the following second-order equation for <p(s): 

(2.5) (n + l)\P 2 (tp)\ n PM = -<p, 

where P k denote linear differential operators (see [33, p. 140]) given by the recursion 
P k +i(<p) = P' k {y) + k)P k { V ), k>0; P (y?) = (p, so that 

p x ((p) = <p' + p 2 &) = (f" + (2/i - iy + n(ji - ly 

P 3 {cp) = ipf" + 3(/z - iy + (3/i 2 - 6/i + 2y + - - 2)<p, 
Pi((p) = <^ (4) + 2(2fi - Z)(ff" + (6/i 2 - 18// + liy 
+ 2(2/i 3 - 9/i 2 + 11/i - 3)<p' + n(n - l)(/i - 2)(n - 3)<p, etc. 

According to f!2.4j) . we are interested in uniformly bounded global solutions ip(s) that 
are well defined as s = lay — > — oo, i.e., as y — > + . The best candidates for such global 
orbits of (12.51) are periodic solutions <^*(s) that are defined for all s6l These describe 
suitable (and, possibly, generic) connections with the interface at s = — oo. The following 
result is proved by shooting as in [TB, § 7.1] and follows the arguments in [231 § 2]. 



Proposition 2.1. For all n > 0, (\2.5\) has a periodic solution of changing sign (p*(s). 
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Periodic oscillatory component, n=0.75 




Periodic oscillatory component, m=2, n=1 




(a) n = 0.75 



(b) n = 1 



Periodic oscillatory component, m=2, n=3 




Periodic oscillatory component, n=5 




(c) n — 3 (d) n = 5 

FIGURE 1 . Convergence to a stable periodic orbit of the ODE (|2.5f) for n 
where ~ 10~ 7 , n = 1, n = 3, and n = 5, with ~ 1CT 2 . 
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There are two open problems: 

(i) uniqueness of the periodic solution <p*(s), and 

(ii) stability y?*(s) as s +oo. 

Numerical evidence answers positively to both questions. Then (i) and (ii) mean a unique 
(up to translation) periodic connection with s = — oo, where the interface is situated. 

The convergence to the unique stable periodic behaviour of (I2.5P is shown in Figure [1] 
for various n = 0.75 (periodic oscillations are of order 1CT 7 ) and n = 5 (order is 10" 2 ). 
Different curves therein correspond to different Cauchy data ^(0), ^p'(0), (p"(0) prescribed 
at y — 0. For n < |, the oscillatory component gets extremely small, so an extra scaling 
is necessary as explained in [T6l § 7.3]. A more accurate passage to the limit n — > in 
the degenerate ODEs such as (12. 5p is presented there in Section 7.6 and in Appendix B. 

Finally, given the periodic <^*( s ) of (12. 5p . as a natural way to approach the interface 
point yo = according to (12.41) . we have that the ODE (12.31) and, asymptotically, (12. 2p . 



Non-periodic oscillatory component, m=2, n=t x1(T s Non-periodic oscillatory component, enlarged, m=2, n=1 




(a) stability of equilibria (b) enlarged behaviour 

FIGURE 2. Non-oscillatory behaviour for the ODE (f!TT|) for n = 1; stability of 
equilibria (|2.8p (a), and enlarged non-periodic behaviour in between, (b). 



admit at the singularity set {/ = 0} 

(2.6) a 2D local asymptotic family with parameters yo and phase shift in s t— > s + sq- 
We also call (12.61) an asymptotic bundle of orbits. 

2.2. Non-oscillatory case A > 0. For A = 1, we have the opposite sign in the ODE 
(2-7) (n + l)\P 2 (<p)\»P s (<p)=<p, 

which admits two constant equilibria 

(2.8) ^ ± = ±[(n + l)(//-2)]^(//-l)]^. 

Figure E^a) shows that as s — ► +oo the equilibria (12.81) are stable (easy to see by lineariza- 
tion). In (b), which gives the enlarged behaviour from (a) close to (p — 0, we observe 
a changing sign orbit, which is not periodic. This behaviour cannot be extended as a 
bounded solution up to the interface at s = — oo. In other similar ODEs, which are in- 
duced by other parabolic PDEs, such behaviour between two equilibria can be periodic; 
cf. [331 P- 143]. 

These results confirm that for A > 0, the TWs are not oscillatory at interfaces, and 
actually such backward propagation via TWs is not possible for almost all (a. a.) initial 
data. More precisely, unlike (I2.6p . for A > 0, the asymptotic family (a bundle) as s — > — oo 
is ID, which is not sufficient for matching purposes (see typical ideas of construction of 
similarity profiles below). 
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3. Blow-up similarity solutions: problem setting and preliminaries 

3.1. ODE reduction. The parabolic PDE (11. 6ft formally possesses the following simi- 
larity solutions describing finite-time blow-up as t — » T~: 

(3.1) u s (x,t) = (T-t)-^f(y), y = x/{T-t)f>, with P = 2 &^l iy 
The rescaled blow-up profile f(y) satisfies the quasilinear fourth-order ODE 

(3.2) A(/) = -(|/T/T-/3y/ / -^T/+i/r 1 / = in R. 
We impose at the origin y = either the symmetry conditions, 

(3.3) f(0) = and /'"(0) = 0, 
or the anti-symmetry ones, 

(3.4) /(0) = and f(0) = 0. 

By a standard local analysis of (13.21) for small / ~ 0, and in view of general results 
on regularity [391 Ch. 1,2] and finite speed of propagation for such degenerate parabolic 
equations [28], a natural setting for the Cauchy problem assumes that, for p e (1, n+ 1], 

(3.5) f(y) is sufficiently smooth and compactly supported. 

The actual regularity of f(y) close to interfaces has been determined in the previous 
section. 

For p > n + 1, the asymptotic analysis shows that the solutions are not compactly 
supported. Note that equation (13. 2p possesses the constant equilibria 



(3.6) ±Up) = ±{p 



p-i 



3.2. Blow-up self-similar profiles: preliminaries. We next study solvability of the 
ODE (13. 2p in R. First of all, the local interface analysis from Section [2] applies to (13. 2\i . 
Indeed, close to the interface point y = y > of the similarity profile f(y), the ODE 
(13. 2p for p < n + 1 contains the same leading terms as in (12. 3p and other linear two are 
negligible as y — > y . 

For p — n + 1, where (3 = 0, the leading terms close to the interface are 

-(l/T/")"-£/ = o. 

This gives solutions (12.4ft with another exponent 

2(n+2) 
• n ' 

and a fourth-order ODE for <p{s), which admits a periodic solution v 9 *( s ); see examples 
in [331 Ch. 3-5]. 

It is important that, taking into account the local result (12.61) and bearing in mind the 
two conditions (13.31) or (13. 4p yield two algebraic equations for two parameters {y , s } of 
the bundle. Therefore, we expect that 

(3.7) there exists not more than a countable set of solutions. 
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Note that this assumes a certain analyticity of the dependence on parameters in the 
degenerate ODE (13.21) . which is not easy to prove. In particular, relative to the parameter 
p > 1, we can expect at most a countable set of p-branches of solutions. This is true for 
the linear case n = and p = 1; see below. 

3.3. Fundamental solution and necessary spectral properties. Here we review 
some properties of differential operators in the linear case n = 0. Consider the linear 
bi-harmonic equation 

(3.8) u t = —u xxxx in Rx R + . 
Its fundamental solution has the form 

(3.9) b{x,t)=r 1 *F(y), y = x/t 1 *, 

where the rescaled kernel F is the unique radial solution of the ODE 

(3.10) BF = -F (4) + \yF' + \F = in M, with J F dy = 1. 

R 

On integration once, we obtain a third-order equation, 

(3.11) —F'" + \yF = in 1. 

The kernel F = F(\y\) is radial, has exponential decay, oscillates as \y\ — > oo, and 

(3.12) \F(y)\ < De"* |4/3 in R, 

for a positive constant D and d = 3 • 2~ n / 3 ; see [HI p. 46]. The necessary spectral 
properties of the linear non self-adjoint operator B and the corresponding adjoint operator 
B* are of importance in the asymptotic analysis and are explained in [13] for general 2mth- 
order operators (see also [H2, § 4]). In particular, B has a discrete (point) spectrum er(B) 
in a weighted space Lp(M), with p{y) = e a '^' 4/3 , a G (0, 2d) is a constant, 

(3.13) a(B) = {Xt = -{, Z = 0,1,2,...}. 
The corresponding eigenf unctions are given by 

(3.14) Mv) = L $r-F®(v), l = 0,1,2,.... 

The adjoint operator 

(3.15) B* = —Dy — | yD y 

has the same spectrum (13.131) and polynomial eigenfunctions 

(3.16) rM = ^ E J JT^V, / = 0,1,2,..., 

i=o 

which form a complete subset in L^»(R), where p* = -. As B, the adjoint operator B* has 

compact resolvent (B* — A/) -1 . It is not difficult to see by integration by parts that the 
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eigenf unctions (13.141) are orthonormal to polynomial eigenfunctions {ipf} of the adjoint 
operator B*, so 

(3.17) (i>i,i>t)=6i k , 

where (-, •) denotes the standard (dual) scalar product in L 2 (IR). 

3.4. Countable set of similarity solutions for n = 0, p = 1. Performing in the 
equation (11.81) the change 

(3.18) u(x,t) = e t w(x,t) 

reduces it to the pure bi-harmonic equation (I3.8P for w(x,t). By the scaling as for the 
fundamental solution b(x,t) in (13.91) . 

(3.19) w(x, t) = t~*v(y 7 r), y — x/t^, r = lnt, 

we obtain the rescaled equation with the B in (13. 1Q|) having eigenfunctions (13.141) . so 

(3.20) v T = Bv 3 vi{y,r) = e x ' T ^{y). 

Setting A; = — | as in (I3.13P and t = e T , we obtain a countable set of different asymptotic 
patterns for the linear PDE (11.81) corresponding to n = and p = 1: 

(3.21) Ul (x,t) = e"* r^^(-^), I = 0, 1,2, ... . 

It turns out that the blow-up similarity patterns (13. ip can be deformed as n — > and 
p — > 1 to those in (13.211) (though entirely rigorous proof is very difficult and not fully 
completed for such degenerate equations, as will happen for some other related homo- 
topy questions). Then (13.211) suggests that there exists a countable number of branches 
{fi(y, n,p)}, which appear from the branching point {n = 0, p = 1} according to classic 
theory, [371 § 56]. We claim that the above two (linear for n — 0, p — 1 and nonlinear for 
n > 0, p > 1) asymptotic problems admit a continuous homotopic connection as n — > 0, 
p — > 1, so that, after necessary scaling, (I3.2ip is obtained in the limit from nonlinear 
eigenfunctions. For such ODEs, this reduces to a matched asymptotic expansion analysis, 
which is rather technical and is not studied here. 

What is key for the future study is that the oscillatory behaviour of linear patterns in 
(13.211) is then inherited by nonlinear blow-up patterns at least for small n > and p > 1. 
This shows once more that similarity profiles f(y) corresponding to the Cauchy problem 
must be oscillatory near interfaces. Homotopy approaches can play a role for specifying 
correct settings of the Cauchy problem for variety of nonlinear PDEs with non-smooth 
or singular coefficients, if they share the same homotopy class with a well-posed linear 
equation; see [HI Ch. 8]. 

It follows from the ODE (J22D that 

(3.22) ||/IU~/.(p) = (p-ir^ T -^+oo as p - 1+, 

so the divergence (in fact, towards the rescaled linear problem) is exponentially fast. 
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4. Regional blow-up profiles for p — n + l: variational approach 



We begin with the special case p — n+l, where (3 = in ( 13. ip (so y = x) and f(y) in 
(13. 2p solves an autonomous fourth-order ODE of the form 

(4.1) A(/) = -(|/T/T-J/ + l/l B / = in R. 

This is a variational problem that can be studied in greater detail. Later on, we apply 
these patterns and classification for p = n + 1 in neighbouring parameter ranges p > n + 1 
and p < n + 1 by using a natural idea of p-branches of solutions. 
For convenience, we perform in (14. ip an extra scaling 

(4.2) /=(I)»F =>. -(\F"\ n F")"-F + \F\ n F = in R. 
For any n > 0, this equation admits three constant equilibria 

F = -l, 0, 1. 

4.1. Variational setting and compactly supported solutions. Operators involved 
in the ODE (14.21) are potential in L 2 , so the problem admits a variational setting and 
solutions can be obtained as critical points of a C l functional of the form 

(4.3) E(F) = -^g / \F'T 2 dy - f / F 2 dy + ^ / |F|«+ 2 dy. 

Then we are looking for critical points in W 2 n+2 (R) n L 2 (R) n L n+2 (R). For compactly 
supported solutions (see below), we choose a sufficiently large interval B R = (-R, R) and 
consider the variational problem for (14. 3 p in W2q 2 (Br), assuming Dirichlet boundary 
conditions at the end points 8Br = {±i?}. By Sobolev embedding theorem, W£q 2 (Bji) 
is compactly embedded into L 2 (Br) and L n+2 {Bfi). Continuity of any bounded solution 
F(y) is guaranteed by Sobolev embedding if 2 (R) C C(R). 

Thus, we will be looking for compactly supported solutions. This demand is associated 
with the well-known fact that the corresponding parabolic flow with the elliptic operator 
as in (Oj) . 

(4.4) w t = -(\w xx \ n w xx ) xx - w + \w\ n w, 

describes processes with finite propagation of interfaces. By energy estimates, such re- 
sults have been proved for a number of quasilinear higher-order parabolic equations with 
potential p-Laplace-type operators; see [H] . Therefore, our blow-up patterns are indeed 
nontrivial compactly supported stationary solutions of (14.41) . Examples of ODE proofs 
via typical energy estimates can be found in [H § 7]. 

Thus, in what follows, to revealing compactly supported patterns F(y), we will pose 
the problem in bounded sufficiently large intervals (-R, R) with Dirichlet data at ±R. 

4.2. L— S theory and direct application of fibering method. The functional (14.31) 
is C 1 , uniformly different iable, and weakly continuous, so we can apply classic Lusternik- 
Schnirel'man (L—S) theory of calculus of variations [3?! § 57] in the form of the fibering 
method [4"4"l 
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According to L-S theory and the fibering approach, the number of critical points of 
the functional (14.31) depends on the category (or genus) of functional subset on which 
the fibering is taking place. The critical points of E(F) are convenient to obtain by the 
spherical fibering in the form 

(4.5) F = r(v)v (r > 0). 

Here r(v) is a scalar functional, and v belongs to a subset in W2o 2 (Br) given by 

(4.6) Ho = {v E W$ 2 (B R ) : H {v) =-J \v"\ n+2 dy + J \v\ n+2 dy = l}. 
Then the new functional 

(4.7) H(r, v) = E{rv) = ^ r ™ +2 - \ r<1 1 v 2 dy 
has the absolute minimum point, where 

H' r = r n+1 - r J v 2 dy = r (v) = (J v 2 dy) ~, 

at which H(r (v), v) = -5^2) r o +2 ( v )- 

Therefore, introducing 

(4.9) H(v) = [- 2 J2±*L H (r (v),v)] ^ = J v 2 dy, 

we arrive at the quadratic, even, non-negative, convex, and uniformly differentiable func- 
tional, to which L-S theory applies, [371 § 57]. Searching for critical points of H in the 
set Ti,o, one needs to estimate the category-genus p of the set Ho- The details on this 
notation and basic results for semilinear equations can be found in Berger p, p. 378]. 
The Morse index q of the quadratic form Q in Theorem 6.7.9 therein is precisely the 
dimension of the space where the corresponding form is negatively definite. This includes 
all the multiplicities of eigenfunctions involved in the corresponding subspace. Note that 
Berger's analysis and most of others are dealing with perturbation theory of linear opera- 
tors, which makes it easier to get the genus of necessary functional sets involved. For the 
quasilinear operators that define the set (14.61) by their potentials, an extra study of genus 
is needed (to be performed below). 

For detecting geometric shapes of patterns, we recall that by the minimax analysis of 
L-S category theory [37J p. 387], [21 p. 368], the critical values {ck} and the corresponding 
critical points {vk} are given by 

(4.10) c k = ini TeMk sup^ H{v), 

where F C TCq are closed sets, and M.k denotes the set of all subsets of the form 

BS k - x C Ho, 

where S k ~ l is a suitable sufficiently smooth {k— l)-dimensional manifold (say, sphere) in 
H,q and B is an odd continuous map. Then each member of Aik is of genus at least k 
(available in TLq). It is also important to remind that the definition of genus [27J P- 385] 
assumes that p{T) = 1, if no component of F U F*, where 

F* = {v: —v G F}, 
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is the reflection of T relative to 0, contains a pair of antipodal points v and v* = —v. 
Furthermore, p{F) — n if each compact subset of JF can be covered by, minimum, n sets 
of genus one. 

According to (14.101) . 

Ci < c 2 < ... < Q , 
where l Q = lo(R) is the category of 7Y satisfying (see below) 

(4.11) l Q (R) -> +oo as R -> oo. 

Roughly speaking, since the dimension of the sets T involved in the construction of Jv[ k 
increases with k, this guarantees that the critical points delivering critical values (14.101) 
are all different. 

4.3. Category of Ho gets arbitrarily large as R — > +oo. It follows from j37J p. 385], 
[21 p. 376] (see also [45]) that according to (14.61) . the category l = p(7i ) of the set Hq 
can be associated with the maximal number K = K(R) of nonlinear eigenvalues A& < 1 
of the corresponding elliptic problem 

(4.12) -(|^ ,/ |V / ) ,/ + A*|Vl ft V' = o, ^eW 2 fi (B R ). 

This problem is solved by L-S theory and gives at least a countable set of critical values 
and different critical points of the positive homogeneous functional 

(4.13) / \v\ n+2 dy on the unit sphere Si = { / \v"\ n+2 dy = l}. 
Indeed, given an eigenfunction ip*. ^ with A& < 1, multiplying (14. 12j) by ipk yields 

- / K| n+2 dy + f \M n+2 ^ = (1 - A fc ) / |^| n+2 dy>0 

^ = 5 fe Vfc e Wo, |5 fc | n+2 = [(1 - A*.) / |V> fc | n+2 dy]~\ 

where 5^ > is the necessary normalization factor. By L-S theory, all such nonlinear 
eigenfunctions are different (since correspond to different critical values of the functional), 
so that all of them {t/^, k = 1, K} are linearly independent. In order to estimate the 
genus of 7^o, we take their linear combination 

(4.14) v = C 1 $ 1 + ... + C K $ K e'H , 

so on substitution into the functional in (14.61) we get the following algebraic equation for 
the coefficients C = {C h C K } G R K : 

(4.15) G(C) = -J \C^'{ + ... + C K ^\ n+2 dy + J |C^ a + ... + C K $ K \ n + 2 dy = 1, 
which is an equation of a surface Ck in M. K being symmetric under the reflection 

(4.16) C i — ► -C. 

One can see that, by construction of the normalized eigenfunctions ipk, for any fixed 
fc = 1,2, 

(4.17) G(C) = \C k \ n+2 (l + o(l)) as C k -> oo. 
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It is not difficult to see (using the variational and extremal nature of nonlinear eigenfunc- 
tions) that Ck contains a simple closed connected component, which, in view of (I4.16p . 
is homotopic to the unit sphere S K ~ 1 in M> K . By the "additivity" properties of the genus, 
this implies that 

(4.18) p(H )>K(R)-l. 

We do not know whether this estimate is sharp: optimal estimates of the category (genus) 
of the sets and even multiplicity of nonlinear eigenf unctions for such functionals compose 
a difficult open problem, which persists even for classic p-Laplacian operators as in (11. ip . 
Since the dependence of the spectrum on the length R for (I4.12p is, by simple scaling, 

(4.19) \ k (R) = R' 4 - 2n \ k (l) -> 0+ as R -> oo, k = 0, 1, 2, ... , 

we have that the category p(7io) can be arbitrarily large for R ^> 1, and (14. lip holds: 
Proposition 4.1. The ODE problem has at least a countable set of different so- 



lutions denoted by {Fi, I > 0}, and each one Fi(y) is obtained as a critical point of the 
functional l \4-3§ in ^ (^) w ^ sufficiently large R = R(l) > 0. 

4.4. First basic pattern and local structure of zeros. Let us present numerical 
results concerning existence and multiplicity of solutions for equation (14.21) . In Figure [3j 
we show the first basic pattern for (14 .2p called the F (y) for various n G [0.1,0.7]. These 
profiles are constructed by MatLab by using a natural regularization in the singular term, 

(4.20) -[(e 2 + {F") 2 )^F"]" - F +\F\ n F = Q in R (e > 0). 

Here, the regularization parameter e and both tolerances in the bvp4c solver, typically, 
take the values 

(4.21) e = 10~ 2 or 10~ 3 and Tols = 10~ 3 or 10" 4 . 

For n > 0.5, convergence gets rather slow. For n < 0.7, the global structure of blow-up 
profiles (excluding their fine zero structure, see below) is stable with respect to reasonable 
variations of e and Tols. In fact, this reflects the structural stability of first basic blow- 
up patterns, which the author observed in dozens of other nonlinear parabolic models 
with blow-up. Note that proving stability even in the linearized setting involves non self- 
adjoint operators with non-constant coefficients that leads to several technical difficulties 
and remains open. On the other hand, for n > 1, i.e., for strongly nonlinear diffusion 
operators in (I4.20p . we did not get reliable enough numerical results with the necessary 
accuracy, so we will avoid using such cases for further illustrations. 

Incidentally, this makes it possible to reveal some features of the local structure of 
multiple zeros close to the interface. Figure H] shows how the zero structure of profiles 
F from Figure repeats itself in a "self-similar manner" from one zero to another in 
the usual linear scale. In Figure (b), a "discrete", piece- wise continuous structure for 
n = 0.5 is already revealed, and this is the best we have been able to achieve numerically. 
However, this makes no problem, since the accuracy 10~ 3 achieved in (b) is already in 
agreement with parameters in (I4.2ip . so further improvements make no practical sense. In 
addition, this shows that the discrete and continuous solutions of this difficult variational 

16 



Basic pattern F Q (y) 




y 

FIGURE 3. The first solution F (y) of (TOJ) for various n. 

problem remain very similar even for the present rough meshes, when the discrete features 
become clearly visually observable (as usual, it is a key fact for such numerics). 

Further revealing zero structure and eventually the behaviour such as (12.41) as s = 
Ho ~ V ~~ > + cannot be reliable done in the parameter range (14.211) . In pT5], [16] . for similar 
thin film models, this demanded e and Tols to achieve at least 10~ 12 , which is not possible 
for the current model in view of slow convergence for higher-order p-Laplacians. It is also 
quite a challenge to detect numerically the free-boundary point. The main difficulty is 
to distinguish the nonlinear oscillations via (12.41) and the linear ones in the "linearized 
area", where (I4.20p implies an exponential behaviour for y ^> 1 governed by the ODE 

(4.22) = -£~ n F + ... => F(y) ~e-^ £ ~ n/i cos (^e-%y + c), 

where c is a constant. Actually, we saw not more than first 1-3 nonlinear zeros of the 
type (12 .4p and the rest of zeros corresponded to the linear behaviour (I4.22p . 

4.5. Basic countable family: approximate Sturm's property. In FigureEJ we show 
the basic family denoted by 

{F h l = 0,1,2,...} 

of solutions of (14. 2 p for n = 0.2. This family is connected with the application of L-S 
and fibering theory; see [2UI [3D] . Each profile Fi(y) has / + 1 "dominant" extrema and 
I "transversal" (not from the tail) zeros; see [221 § 5], [3DI § 5], and [23 § 4] for further 
details. It is important that 

all the internal zeros of Fi(y) are transversal, 
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Basic pattern F Q (y); enlarged zero structure _ 4 Basic pattern F Q (y): enlarged zero structure 
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FIGURE 4. Enlarged zero structure of the profile Fo(y) in the linear scale. 



excluding the oscillatory end points of the support. In other words, each profile F) is 
approximately obtained by a simple "interaction" (gluing together) of I + 1 copies of 
the first pattern ±F taking with necessary signs. Such a gluing of oscillatory tails is 
illustrated in Figures [6] and U\ see also further comments below. There is some analytic 
evidence [30, § 5] that exactly this basic family {Fi} is obtained by the classic Lusternik- 
Schnirel'man construction of critical points of the reduced functional ( 14. 9p . A rigorous 
definition of gluing assumes formation of all the internal transversal zeros while the outer 
ones at the end point of the support are the only ones that remain oscillatory according 
to the behaviour (12. 4p with the periodic orbit tp = ip*(ln(y — y)). Some question of global 
behaviour of such patterns Fi(y) for large / remain open, [291 130] . 

Let us forget for a moment about the complicated oscillatory structure of solutions 
near interfaces, where an infinite number of extrema and zeros occur. Then the dominant 
geometry of profiles in Figure [5] looks like it approximately obeys Sturm's classic zero set 
property, which is true rigorously for the case m = 1 only, i.e., for the second-order ODE 

n 

(4.23) F" = -F+\F\ n+1 F in R. 

For (I4.23p . the basic family {FJ} is explicitly constructed by direct gluing together simple 
patterns ±F given explicitly; see [331 P- 168]. Therefore, each Fi consists of precisely l + l 
patterns (with signs ±Fo), so that Sturm's property is clearly true by direct application 
of L-S category theory. 

4.6. Countable family of {F , F }-gluing. Further patterns to be introduced do not 
exhibit as clear a "dominated" Sturm property and are associated with a double fibering 
technique where both the Cartesian and spherical representations of critical points are 
used; see [291 § 3], [30J § 3]. Let us present some explanations. 

The nonlinear interaction of the two first patterns Fo(y) leads to a new family of 
profiles. In Figure [6] for n = 0.2, we show the first six profiles from this family denoted 
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The first basic profile F (y); n=0.2 




The second (dipole) basic profile F^y); n=0.2 




(a) F (y) 

The third basic profile F (y); n=0.2 
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(c) F 2 (y) 

The fifth basic profile F (y); n=0.2 




(e) F 4 (y) 



(b) F 1 (y) 

The fourth basic profile F (y); n=0.2 




(d) F 3 (y) 



The sixth basic profile F (y); n=0.2 
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(f) F 5 (y) 



FIGURE 5. The first six patterns of the basic family {Fi} of the ODE (fl~2l) for 
n = 0.2. 
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by {F + 2 t k,+2}- I n the refined zero structure of the last profile in (b), we already see some 
numerical effects of rather rough meshes, which again do not deny the sufficient overall 
quality of numerics. In each function F + 2,k,+2 the multiindex 

a = {+2,/c,+2}, 

from left to right denotes: +2 means two intersections with the equilibrium +1, then next 
k intersections with zero, and final +2 stands again for two intersections with +1. Later 
on, we will use such a multiindex notation to classify other patterns obtained. 

As a general rule, we point out again that any finite gluing of a pair of patterns ±Fo(y) 
actually means that all internal zeros become transversal. Note that this and all the 
Figures involved are not enough to explain the essence of this complicated and math- 
ematically not fully understood procedure for non-homotopic variational problems |31j . 
The resulting patterns have zeros of infinite order only at the end points of its support. 

In view of the infinite oscillatory character of Fo(y) at the interfaces, we expect that 
the family {F +2t k,+2} is countable, and such functions exist for any even k = 0,2,4, ... . 
Then k = +oo corresponds to the non-interacting pair 

(4.24) F (y + yo) + F (y - y ), where supp F (y) = [-y , y ]. 

It is expected that there exist various triple {F , F , F } and any multiple interactions 
{F , F } of k single profiles, with different distributions of zeros between any pair of 
neighbours (proof is an open problem). 

4.7. Countable family of {— F , F }-gluing. We now describe the interaction of —F (y) 
with Fo(y). In Figure [7] for the case n = 0.2 (which is convenient in terms of rather fast 
convergence of the numerical method employed), we show the first profiles from this family 
denoted by {F_ 2i a;,+2}, where for the multiindex a = {— 2, +2}, the first number —2 
means two intersections with the equilibrium —1, etc. It can be seen that the first two 
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First profiles from family F_ 2k+2 ; n=0.2 
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FIGURE 7. First four patterns from the family {F_2,k,+2} of the {— Fo,Fo}- 
gluing; n = 0.2. 

profiles belong to the same class .F_2,i,2j i>e., both have a single zero for ?/ w 0. The 
last solution shown is F_ 2)5j+2 . Again, we expect that the family {F_ 2) fc )+2 } is countable, 
and such functions exist for any odd k = 1,3,5,..., and k = +oo corresponds to the 
non-interacting pair 

(4.25) -F (y + y ) + F (y -y ). 

We expect that there exist families of an arbitrary number of gluing {±Fo, ±-?o, ±Fo} 
consisting of any k > 2 members (again an open problem). 

4.8. Periodic solutions in R as new types of oscillations about ±1. Before in- 
troducing new types of patterns, we need to describe other non-compactly supported 
solutions in R. As a variational problem, equation (14. 2p admits infinitely many periodic 
solutions; see e.g., [101 Ch. 8]. Figure [8] for n = 0.2 reveals unstable periodic solutions 
obtained by shooting from the origin with various Cauchy data at y — 0. In (b), the pe- 
riodic orbit F*(y) is oscillating about the equilibrium F = — 1. It turns out that precisely 
the periodic orbit F*(y) in (a) with the range 

(4.26) min F,(y) = 0.4135..., max F m (y) = 1.4085... (n = 0.2) 

plays an important part in the construction of other families of compactly supported 
patterns. Namely, all the variety of solutions of (14. 2 p that have oscillations about equilibria 
±1 are close to ±F*(y) there. 

4.9. Family {F + 2k}- Such functions F + 2k for k > 1 have 2k intersections with the single 
equilibrium +1 only and have a clear "almost" periodic structure of oscillations about. 
The number of intersections denoted by +2k gives an extra Strum index to such a pattern. 
In this notation, for k = 1, we have 

F+2 = F . 
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Convergence to periodic solution about equilibrium F=1 ; n=0.2 Convergence to periodic solution about equilibrium F=-1 ; n=0.2 




(a) about F = 1: F(0) = 0.4135 (b) about F = -1: F(Q) = 1.4085 

FIGURE 8. Examples of convergence to periodic solutions of the ODE (|4.2I) for 
n = 0.2; about F = 1 (a) and F = -1 (b). 



Profile F (y) and F (y); n=0.2 
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FIGURE 9. The first profiles F +4 (y) and F +6 (y) from the family {F +2 k, k>2}; 
n = 0.2. 

Two profiles and F + q are shown in Figure [9] for n = 0.2. The further profile 
^+4,1,-2,1,4-4(2/) comprising two sub-structures F +4 from the family {F +2 fe} is shown in 
Figure [10] by the boldface line. 

4.10. More complicated patterns: towards chaotic structures. By combining the 
above rather simple families of patterns, we claim that a pattern (more precisely, a class 
of patterns) with an arbitrary admissible multiindex of any length 



(4.27) 



a = {±01, 02, ±03, &4, •••,±07} 
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Pro,ile F + 4,i,-2,i, + 4<y> ; n=0 - 2 




y 



FIGURE 10. The profile F+i 1-21+4(2/) (the boldface line) from the family 
F + f tj i j -m ! i : +k of solutions of (|4.2p ; n = 0.2. For comparison, the profiles F + ^{y) 
and F + e(y) from Figure [9] are presented. 

can be constructed. For example, in Figure [TU a single complicated pattern with 
(4.28) o = {-2, 2, -2, 1 + 2, 2, +2, 1, -8, 1, +2} 

is given. The computation (the convergence is rather slow for this type of p-Laplacian 
operators) is performed for n = 0.2 as usual. We must admit that, as it is seen from Figure 
[TTl the iterations have not been properly converged within the parameter range (I4.2ip . 
but we can guarantee the convergence with the accuracy at least ~ 10 _1 . This is not that 
bad, since such patterns are not structurally stable and have multi-dimensional unstable 
manifolds, so the convergence must be extremely slow. It is worth mentioning that, using 
bvp4c solver, this computation took a few hours with the maximal number of 75000 points 
on the interval (—50,50). Nevertheless, regardless such a lack of accuracy, we are sure 
that such complicated critical point profiles really exist, since we have seen a lot of those 
in other similar (and simpler numerically) higher-order variational problems [29, 30J that 
were not associated with such awkward and strongly degenerate operator as p-Laplacian 
ones. Special more conservative and divergent numerical techniques are necessary for 
tackling higher-order p-Laplacian operators in the ODEs, but here we demonstrate what 
an average (numerically, non-professional) PDE user can extract from standard MatLab 
codes. Theoretically, via the L-S/fibering theory, all those patterns are well defined. 

We claim that the multiindex (14.271) can be rather arbitrary taking finite parts of any 
non-periodic "fraction". Actually, this means chaotic features of the whole family of 
solutions {F a }. In fact, there is no any exiting news in such a chaotic proclamation: 
one can see that even the basic simple countable family {Fi} is indeed chaotic, since the 
choice of the sequence of elementary profiles ±Fq in Fi for I ^> 1 can be arbitrary long 
with an arbitrary sequence {±} of sign changes, thus exhibiting no finite periodic order in 
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FIGURE 11. A complicated pattern F a (y) of the index (|3~28j) for ([OP for n = 0.2. 

the index a. These chaotic types of behaviour are known for other simpler fourth- order 
ODEs with coercive operators and definite homotopic features, [32, p. 198]. 

A variety of complicated patterns of these types for different variational problems as- 
sociated with the PME-type nonlinearities in (11.121) can be found in [29[ § 5.6], [30j § 5.6] 
(see also § 6 in [30J therein for sixth- and eighth-order models) and in [23j § 5]. The 
convergence of standard numerical methods for these problems are much faster, with s, 
Tols up to ~ irr 10 . 

5. Single point blow-up for p > n + 1: P- and Q-type profiles 

We now return to the similarity ODE (13. 2p in the case p > n + 1, which, in view of 
the spatial rescaled variable y in (13.11) . corresponds to single point blow-up. It is key that 
(13. 2p for p 7^ n + 1 is not variational. Formally, solutions of (13. 2p can be traced out by 
shooting and matching procedures, which are too complicated. Instead, we will use a 
continuation in parameters approach, which allows us to predict solutions by using those 
in the variational case p — n + 1. 

5.1. Asymptotics at infinity and single point blow-up. We begin with simpler 
asymptotics of the solutions of (13.21) as y — ► +oo. Unlike the previous case of regional 
blow-up for p = n + 1, for p > n + 1, equation (I3.2p admits non-compactly supported 
solutions with the following behaviour: 



The first term in (15. ip represents an "analytic" part (can be truly analytic for some 
parameters) of the expansion, while the second one gives the essentially "non-analytic" 
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c - n ] 3 > 0. 



part. Such a structure in (15 .ip is usual for saddle-node-type equilibria [HI p. 311], but, 
f or the fourth-order ODE (13. 2p . this expansion does not admit a simple phase-plane 
interpretation. However, existence of such asymptotic expansions can be justified by 
fixed point arguments, which becomes quite a technical issue and is not done here. 

One can see passing to the limit t — > T~ in (13. ip that the first term in the asymptotic 
expansion (15. ip gives the following final-time profile of this single point blow-up for even 
patterns / = f(\y\): 

2(n + 2) 

(5.2) u s {x,T') = C \x\'p-("+ 1 i < oo for all i^O. 
Returning to the asymptotic expansion, we conclude that (15.11) represents 

(5.3) a 2D asymptotic family (bundle) of solutions. 

Hence, the family (15.11) is well suitable for matching with also two symmetry conditions 
at the origin (13.31) . so we expect not more than a countable set of solutions. For first 
patterns, we keep the same notation fi(y) as in Section H] for p = n + 1. 

5.2. Oscillatory behaviour about constant equilibrium In order to predict the 
multiplicity of solutions of (13.21) . we need to study its oscillatory properties. To this end, 
we perform the linearization about the constant equilibrium /* in (13. 6p of the ODE (13. 2p , 

(5.4) f = f* + Y 

formally assuming that \Y\ <C 1 on some bounded intervals. This yields the "linearized" 
nonlinear equation 

(5.5) B n (Y) = -(\YrY")"-{3Y'y + Y = ( (5 = ) ■ 

We are going to study oscillatory, sign-changing properties of solutions of (15.51) for 
various n > 0. Notice that the "linearized" ODE (15.51) remains a difficult fourth-order 
equation. Indeed, in view of the invariance with respect to the group of scalings 

2Q+2) 

Y I— > e ~Y, yv-^ey (e > 0) 

the transformation 

2(n+2) 

Y(y) = y » ^(s), s = \ny, 

reduces (15.51) to an autonomous fourth-order ODE. Setting P(^p) = <p' yields a third-order 
ODE, but further reductions are impossible. Thus, (15.51) cannot be studied on the phase- 
plane in principle; cf. |7J, where, for (11.31) . oscillatory analysis on the phase-plane is a 
convenient and exhaustive tool. 

Therefore, we will need another further investigation of (15. 5p . and we begin with the 
following useful comment: 

Linear case n = 0. Then the quasilinear operator B n in (15.51) becomes linear, 

(5.6) B F = 0, where B = -D\ - \ yD y + / = B* + / (/3 = f) 
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and B* is the adjoint linear operator (|3.15|) . According to its point spectrum (I3.13p . 
equation (15.51) for n = has a non-oscillatory solution, being the eigenfunction ipf(y) for 
I = 4, i.e., 

(5-7) Y(y) = r i (y) = ^(y i + 2A) (n = 0), 

which is an example of a non- oscillatory solution. Nevertheless, one can see from the 
operator (15.61) that the linear ODE (15.51) for n = has other oscillatory solutions with an 
increasing envelope as y — > +oo; see below. 

Quasilinear case n > 0. In Figure [T27a)-(d) . we show that, for any p > n+1, the ODE 
(15. 5p admits infinitely oscillatory solutions with increasing amplitude of oscillations. 

Here, (a) shows linear increasing amplitude of oscillations for n = 0. It is curious 
that such behaviour persists in the nonlinear range n > 0, p > n + 1, so that n = is 
a branching point for (I5.5P from solutions of the linear equation (I5.6p . In (b), we show 
the bounded periodic solution for the variational case p = n + 1 (cf. Figure E]), which 
generates a p-branch of non-periodic patterns for p > n + 1; see (c). This suggests that 
basic blow-up similarity patterns {fi(y)} are expected to exist for p > n + 1 sufficiently 
close to n + 1. Figure fT2T d) shows that for p < n + 1, the amplitude of oscillations becomes 
decreasing, so we expect a single P-type profile fo(y) for p £ (l,n + 1). 

5.3. Basic patterns: first numerical conclusions. For convenience, similar to the 
change in (14.21) . we perform the following scaling in (13.21) for p ^ n + 1: 

(5.8) / = AF, y i — > ay, where A = /* = {p — 1) _ ~, a — (p — 1)^. 
Then F = F(y) solves the equation 

(5.9) -{\F"\ n F")" -(3F'y-F+\F\P- l F = 0, with $ = {p - l)/3 = ^gj 1 , 

which has the scaled equilibria F* = ±1 that are convenient for numerical experiments 
for small n > and p close to 1 + . 

In Figure [T3"| we present the first pattern F (y) for n = 0.2 for p = n + 1 = 1.2 (the 
dotted line for comparison), 1.4, 1.8, 2.2, 2.6, 3, 4, and 6. In particular, it is clearly seen 
that, for larger p, the profiles approach the positive asymptotic behaviour (15.11) for y ^> 1 
with Co > 0, and become strictly positive in KL 

Figure [H] shows the first F (y), the third F 2 {y), and the fifth F^y) P-type patterns for 
n = 0.2 and p = 1.5. In Figure [T5| we demonstrate two profiles Fo(y) and F 2 {y) for the 
same n = 0.2 and p = 2.6. In the last case, the next even profile F±(y) was not detected 
numerically, and this nonexistence will be confirmed later by the p-branching approach. 

5.4. Q-type profiles. It follows from (13. 2p that an oscillatory expansion can be started 
at any finite point y = yo > 0, at which 

(5.10) f(y ) = U f(y) = h for y < y . 

This yields the so-called Q-type solutions; see classification in [271 E] • Then setting 

y^yo + y, y > o, 
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Linear exponential oscillations for n=0 Periodic solution of the "linearized" equation for n=0.2, p=n+1 =1 .2 
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(a) n = (b) p = n + 1 : n = 0.2 



Oscillations for n=0.2, p=1.4 Decreasing oscilllatory solution for p<n+1 : n=0.2, p=1.1 




5 10 15 20 25 5 10 15 



y y 

(c) p > n + 1 : p = 1.4, n = 0.2 (d) p < n + 1 : p = 1.1, n = 0.2 

FIGURE 12. Examples of oscillatory patterns of the "linearized" ODE (|5.5|) . 

we again arrive asymptotically at a linearized equation similar to (15. 5p . 

(5.11) B n (Y) = -{\Y"\ n Y")" - \ Y' = 0, A = [3y > 0, 

so on integration we obtain the TW equation (12.31) . where the constant Ao is scaled out. 
Therefore, we use the change (12 ,4p to get the oscillatory ODE (12.51) . According to (I2.6p . 
this gives a 2D asymptotic family to be matched with the bundle ( 15.11) at infinity. 

Analytically, as well as numerically, the problem of existence of a countable subset of 
such Q-type similarity profiles is more difficult. Figure [TBI shows the first Q-type profile 
Fq(h) for n = 0.2 for p = 1.5. The convergence here is slower and we do not succeed 
in getting other Q profiles. For the sake of comparison, we also present here P solutions 
F (y) and F 2 (y). 



5.5. On branching of solutions from variational critical points. Consider the ODE 
(15. 9p from the point of view of a perturbation approach. For p = n + 1, i.e., for (3 = 0, 
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Basic profile F (y) for n=0.2, p e [1 .2,6] 




FIGURE 13. The first basic single point blow-up patterns Fo(y) of (|5,9p for 
n = 0.2 and various p € [1.2, 6]. 



P-type profiles F F 2 , F 4 for n=0.2, p= 1 .5 




FIGURE 14. Three single point blow-up patterns i*b(y)> ^2(2/)) and -£4(2/) of 
n = 0.2, p = 1.5. 



the ODE has been studied in SectionHl Setting e = p— (n+ 1) and assuming that |e| > 
is sufficiently small, we write (I3.2p in the form 



(5.12) F(/) = -(|/T/T-i / + !/!"/ 



2(n+2)(p-l) 



yf-^f + \f\ n f(i-\fn 



2.x 



P-type profiles F Q and F 2 for n=0.2, p= 2.6 
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FIGURE 15. Two patterns F (y) and F 2 (y) of (^Mj; n = 0.2, p = 2.6. 



Q-type profile Fr'(y) for n=0.2, p= 1 .5 




y 



FIGURE 16. The first Q-type solution Ftf (y) (the boldface line) of the ODE 
(H2H); n = 0.2, p = 1.5. 

On the right-hand side, the key perturbation term satisfies 

(5.13) 0(/,e) = |/rY(l-|/n-O as e -> 

and is at least continuous at the singular point / = 0, e = 0. On the left-hand side of 
(15.121) . we have the variational operator from Section HI Thus, the non-perturbed problem 
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(15.121) for e = admits families of solutions described in Section 0] (see also extra details 
on other families of solutions in [2H § 5] and [301 § 5,6]). Therefore, classic perturbation 
and branching theory [TQl ED, E7j suggests (and does not prove since g ^ C 1 at (0,0)) 
that, under natural hypotheses, the variational problem for e = generates a countable 
family of p-branches, which can be extended for some sufficiently small |e| > 0. 

The analysis of bifurcation, branching, and continuous extensions should be performed 
for the equivalent to fl 5 . 1 2 [) integral equation with Hammerstein compact operators; see 
typical examples in [8j [251 [36], where similar perturbation problems for blow-up and 
global patterns were investigated. Namely, denoting by D = (Dy)^ 1 the inverse of Dy in 
a sufficiently large interval (-R, R) (for p < n + 1) and in K (for p > n + 1), we write 
(I5TT2D as follows: 

— n 

f = A(f, e) = D D/i " +1 D/il, where 

(5.14) K ' U 1 J 

Kf, e) = y/' - ^ f + " l/i £ ) - ( - J / + l/l"/)- 

For e = 0, this gives the integral equation with a potential operator, which is equivalent 
to the differential one studied before. 

As customary, parameter p-branches of solutions of integral equations with compact 
operators are fully extensible and can end up either at a singularity point or at another 
bifurcation values; see [10] for a modern sounding of such results. An efficient way to 
prove branching is using degree-index theory, which establishes branching from an isolated 
solution^, say, the first one fo for simplicity, from the branching point e = provided its 
index satisfies [371 P- 353] 

(5.15) 7^0. 

For differentiable operators, where the index ind(/ , 1 — A'(/o, 0)) is equal to the rotation 
of the vector field / — A'(/ , 0), there are special techniques for its calculations; see [37J 
Ch. 20,24]. 

However, in view of non-potential structure in (15.141) . and of a non-sufficient (in the 
usual sense) differentiability of (15 . 1 3[) at (0,0), it is convenient to use other alternatives 
of bifurcation-branching theory without direct differentiability hypotheses that assume 
sufficient regularity of the perturbations; cf. [121 § 28]. Namely, in view of our difficulties 
with the differentiability, using Theorem 28.1 in [121 P- 381] replaces differentiability by 
a slightly weaker control of smallness of nonlinear terms in a neighbourhood of (0, 0). As 
usual, the key principle of branching is that it occurs if the corresponding eigenvalue has 
odd multiplicity (the even multiplicity case needs an additional treatment, which is also 
a routine procedure not to be treated here); see further comments below. Justification 
of branching phenomena for such quasilinear degenerate p-Laplacian operators (includ- 
ing also questions of compactness) in the present problem needs further deeper analysis 
and more involved functional topology/constraints. Therefore, most of further analytical 
conclusions remain formal and are open problems. Nevertheless, it turns out and will be 



Proving that a given solution is isolated is also a difficult problem, especially for p-Laplacian operators. 
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checked numerically, the predicted branching behaviour from p = n + 1 actually occurs, 
and thus becomes a key tool of the proposed study of non- variational problems at hand. 

Thus, anyway, for actual applications, one needs to know the spectrum of the self- 
adjoint operator F'(/o). We have that ( 15. 14ft contains nonlinearities that are hardly dif- 
ferentiable at 0, so these applications lead to difficult technical problems, where branching 
analysis from oscillatory profiles fo(y) with interior transversal zeros needs taking into ac- 
count more complicated "functional topology" . Anyway, continuing the application, it is 
worth mentioning that 

A = (with the eigenfunction ip ~ /q) 

is always an eigenvalue of F'(/o). Actually, this corresponds to the invariance of the 
original PDE (I1.6P relative to the group of translations with the infinitesimal generator 
D y . In particular, this implies that A = 1 is an eigenvalue of A'(/o, 0), so this corresponds 
to the critical case, where computing of the index is more difficult and is performed as in 
[23 § 24]. It is more important that 

A = 1 (with the eigenfunction tpi ~ / ) 

is also an eigenvalue of F'(/ ) (this is associated with the group of translations with the 
generator D t ). In other words, the index condition (15.151) or smallness assumption via 
[T2l Th. 28.1] need special additional treatment associated with spectral properties of the 
linearized operator F'(/o). In this connection, we conjecture that branching at p = n + 1 
is valid in appropriate functional setting, and there exist continuous e-curves, from any 
profiles from the family of basic patterns {//, / > 0} constructed in Section I4~5l 

Finally, let us note another convenient (but not that efficient) way to use Schauder's 
Theorem applied to (15. 14j) to get solutions of (15.121) for small e > and to trace out 
p-branches of the suitable profiles. This approach effectively applies in the case of porous 
medium operators as in f 1 1 . 1 2 [) . [23, § 6]. 

5.6. Numerical construction of p-branches. We recall that (I5.9P is not variational, 
and we are going to use a certain continuity feature concerning the limit p — > n + 1. 

The basic p -branch of F (y) of the simplest shape (as well as the pi-one for I = 1) exists 
for all p > 1. In Figure [TTJ, we show the first p-branch of F in (a) and the deformation of 
the profiles Fo(y) in (b), for n = 0.2 and p G (1.05, 6.15). This branch is extended to the 
global blow-up case p < n+1, to be discussed next in Section [61 We expect that, as usual 
in blow-up analysis, this first po branch is composed from structurally stable solutions 
and hence represents the generic blow-up behaviour for the parabolic PDE (11.61) . 

Deformation with p of Fi(y) (the part for y > is shown only) for p slightly above 
the variational exponent p — 1 +n — 1.2 for n = 0.2 is presented in Figure O Further 
extension of this branch beyond p = 1.218 leads to strong numerical instabilities that 
possibly reflects the actual nonexistence of such solutions far away from p — n+1. 

Concerning other, more complicated profiles for p = n + 1, such as F + i(y) and others 
containing structures shown in Figures M, El and[10l numerical results suggest that some 
of them cannot be extended for p > n + 1 (then any version of (I5.15P is not valid). 
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Nevertheless, for F + 4 this is not the case; see Figure [THl where the extension is shown to 
exist for all p > n + 1 and that 

(5.16) 1 1 -^+4 1 1 oo — ^ 1 + as p — > +oo. 

It seems that all the global p-branches satisfy (15.161) ; cf. an analogous result in [25] for 
global similarity solutions. We expect that a similar p-branch of _F +4 is originated at a 
saddle-node bifurcation for some G (l,n + 1), at which it appears together with the 
p-branch of the profiles -F+2,2+2; see further comments below. 
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p-branch of F +4 (y) for p e [1 .2, 5], n=0.2 



Deformation of F (y) for p e (1 .2,5), n=0.2 




The p-branches can connect various profiles, with rather obscure understanding of pos- 
sible geometry of such branches and their saddle-node bifurcation (turning) points. For 
p = n + the questions on connections with respect to regularization parameters as 
in (I4.20p are addressed in [31] posing problems of homotopy classification of patterns in 
variational problems and approximate "Sturm's Index" of solutions. 

6. On global blow-up similarity profiles for p e (l,n+ 1) 

6.1. Local oscillatory behaviour close to interfaces remains the same. Indeed, 
the ODE (HL2D now reads for / « as 

-(\f"\ n f")"-Pyf' + - = o (p<o), 

and reflecting near interface y i— > y — y , on integration for small y > 0, we have 

(\f"\ n f")' = Pyof + .... 

This is precisely (12.31) with A = /% < to be reduced to —1 by scaling. Hence, for 
p G (l,7i + l), the similarity profiles are equally oscillatory near interfaces as for p = 
n + 1. Therefore, according to (12. 6p . this local 2D asymptotic family looks sufficient to be 
matched with two symmetry boundary conditions (13.31) (or (13. 4p ) at the origin, though 
the proof of existence remains open. 

6.2. On similarity profiles and p-branches. For 1 < p < n+1, the rescaled ODE (15.91) 
is more difficult to solve numerically than for p > n + 1. Figure [201 shows deformation of 
similarity profiles Fo(y) for n = 0.2 and p e [1.05, 1.2]. We observe an easy visible growth 
of solutions as p — > 1~. Structurally, the first basic profile i*b(y) remains of a similar 
geometric shape as in the variational case p = n + 1 = 1.2. 
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First basic profile F n (y) for n=0.2 and p < n+1 = 1 .2 
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-10 -5 5 10 

y 

FIGURE 20. Deformation of the first basic profile F (y) of the ODE ([BU]) for 
n = 0.2 and p = 1.2, 1.15, 1.1, 1.05. 

A part of the corresponding p-branch of profiles Fq was shown earlier in Figure fTWa). 
More detailed and sharp results are presented in Figure |2~T1 for p e [1.023, 1.2], n = 0.2, 
where we used branching from the variational profile for p = 1.2 (with the step size 
Ap = — 10 -3 ). From (a), we definitely observe that this pn-branch is going to blow- 
up as p — > 1~, as suggested before. Note that the po-branch is expected to consist of 
asymptotically (structurally) stable blow-up profiles F (y), but we cannot prove this even 
in the linearized approximation. The linearized operator is a difficult non-self-adjoint one 
with unknown spectrum and proper functional setting. 

The next p-branch of dipole-like profiles Fi(y) is shown in Figure l2"2Ta). together with 
the deformation (b) of the functions Fy{y). It is seen that this p-branch is global and 
blows up as p —* 1~. 

We claim that the p-branches of the basic similarity profiles {fi{y)} (q-v. (I5.8P ) are 
extended up to p = 1~, with a blow-up behaviour as in (13.221) . We expect that these 
p-branches can be connected as n — > with those predicted by the linear problem with 
patterns (13.211) . We refer to [361 P- 1090] for an example of such an analysis. For instance, 
in Figure [231 we present the p-branch and the corresponding deformation of the third 
basic profile F 2 (y), which for p = n + 1 is given in Figure E^c) (with the opposite sign). 

Also, a principal fact of existence of the p-branch of the non-basic profiles is explained 
in Figure [2H where a local p-branch of F +i (y) (see Figure M for p — n + 1) is shown to 
exist for p < n + 1 = 1.2 for n = 0.2. 

It is key that this branch cannot be extended for all 1 < p < n + 1. We expect that, 
as p < n + 1 decreases, the p-branch of F + ^ meets the p-branch of the "geometrically 
similar" profile F+2,2,+2 shown in Figure [6](a) (both have two dominant maxima and a 
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single minimum in between) in a turning saddle-node bifurcation point (another branch 
scenario is also possible, [231 §7]). Such scenarios were detected in variational problems; 
see [31]. In the present non- variational case, both analytic and even a reliable numerical 
description of such bifurcations become much more difficult and still obscure. 

Acknowledgements. The author would like to thank S.I. Pohozaev for a discussion 
on genus theory. 
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Deformation of F +4 (y) for p e [1.185,1.2], n=0.2 

1.5 i 1 1 1 1 1 —;— 




y 

FIGURE 24. The deformation of F +A {y) of the ODE ([53]) for n = 0.2 and 
p = 1.2, 1.195, 1.19, 1.185. 
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